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Calculations of the Casimir energy for spherical geometries which are based on integrations of 
the stress tensor are critically examined. It is shown that despite their apparent agreement with 
numerical results obtained from mode summation methods, they contain a number of serious errors. 
Specifically, these include (1) an improper application of the stress tensor to spherical boundaries, (2) 
the neglect of pole terms in contour integrations, and (3) the imposition of inappropriate boundary 
conditions upon the relevant propagators. A calculation which is based on the stress tensor and which 
- - - avoids such problems is shown to be possible. It is, however, equivalent to the mode summation 

j»j , method and does not therefore constitute an independent calculation of the Casimir energy. 
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J>-^, In 1948 Casimir [1] first predicted that two infinite parallel plates in vacuum would attract each other. This 

Ci • remarkable result has its origin in the zero point energy of the electromagnetic field. While the latter is highly 

^> ' divergent, the change associated with this quantity for specific plate configurations is found to be finite and thus 

in principle observable. Early work to detect this small effect [2] was characterized by relatively large experimental 

^^ [ uncertainties which left the issue in some doubt. More recent efforts [3] have provided quite remarkable data, but are 

^ . based on a different geometry from that of Casimir. Since a rigorous theoretical calculation has never been carried out 

*vj ' for the latter configuration, there remains room for skepticism as to whether the Casimir effect is as well established 

^ , as is frequently asserted. 

f^ — The extension of Casimir 's result to problems with nonplanar boundaries has been of considerable interest and 

l/~) , fraught with difficulties. It was achieved for the case of the conducting spherical shell by Boyer[4] in a remarkable 

^^ ' but intricate calculation. His result was subsequently verified by various methods [5-7] , including in particular that of 

^^ [ direct mode summation [8-9] . One of the methods developed in ref . 7 consisted of integrating the radial component of 

^m^ • the stress tensor over the bounding surface, a technique which was subsequently applied to the case of the cylinder [10] , 

^y. ' the circle [11], and the Dirichlct problem of a D-dimensional sphere [12]. In view of the increasingly wide application 

" •^ , of this technique, it is clearly of interest to ascertain its validity. 

'~i ■ Since the electromagnetic sphere is the only known case in which a finite Casimir energy can be obtained for 

I , a nonplanar geometry using conventional (e.g., exponential) regularization techniques, it is convenient to use this 

!jT^" case as a specific framework for the present work. One begins by defining the usual stress tensor for the uncoupled 

(— I I electromagnetic field 



X 



^. where g'^'^ = (1,1,1,-1) and 



j^t^u ^ ^A""^^ _ lgM^J^"/3^^^ (1) 



pt.^ = d^A" - ^•'A^' 



with A^ being the usual vector potential [13]. The latter is taken most conveniently to be in the radiation gauge 
dkA'^ — Q,A'^ — 0. The vacuum expectation value of T^'^ can clearly be obtained from appropriate derivatives of the 
time ordered product 

G'^-(x,i;x',0 = *(0|(A'(x,i)A^(x',O)+|0) (2) 

where the latter is to be evaluated in the presence of spherical conducting boundaries at r = a and r ^ R where 
r = |x|. The limit i? — > oo is to be taken as the final step of the calculation. 
The differential equation for C?*-' 

-^' + ^) G''(x,t;x',0 = 5'^(x-x')5(i-i') 

with (5*''(x) the usual transverse delta function is to be solved for < ?■, r' < a and for a < r,r' < R subject to the 
usual causal boundary conditions, namely positive (negative) frequencies for t ~ t' > (t — t' < 0). This is an issue 



of some importance since all previous stress tensor calculations of the Casimir energy for nonplanar boundaries have 
imposed outgoing wave conditions in the region exterior to the surface of interest. While causal boundary conditions 
follow directly from the existence of the vacuum as the state of lowest energy, there is no way to infer outgoing wave 
conditions [14] for the propagator (2). 

The Fourier transform of the propagator is thus inferred from the boundary conditions to have the form 

G^-'-(x,x';^) = V4^^^14-^ (3) 

n " 

where the sum is to be taken over all normalized eigcnfunctions A\^ which satisfy the wave equation 

(V2+c.2)A„(r)=0 

subject to appropriate boundary conditions at r = 0, a, i? together with the transversality condition dkA^ = 0. From 
r-B = rxE = for r = a,R where E^ = F^^ and Bk — ^^kimF''™', one readily obtains the form of these functions. 

Using finiujlJr) and gi„{io\Jr) respectively to denote the TE and TM modes the eigenfunctions are seen to be of the 
form 



and 



A|l(x)-/,„(a>«r)X;„,(0,0) 

where 'Kirn{0,(t>) with / > 1 denotes the vector spherical harmonics [15], and /;„ and gin are those real linear combi- 
nations of the spherical Bcssel functions ji and ni which satisfy the boundary conditions. 

Denoting the interior and exterior radial functions by {fi^,9in) and {fi^,9i^) respectively, it follows that for < 
f ^ CL, fj^ and g^^ are proportional to ji and required to satisfy the conditions 

from which the interior eigenfrequencies W;^ are determined. For a < r < R the relevant radial functions are 
and 

dlni^ln n = ClnJli^ln ^ + dlnni[uJi^ r) 



where aim h-m cin, and din are constants. The ratios bin/ain and din/cin together with the exterior eigenfrequencies 

dr 



'^in (•^ — li 2) are fixed by the vanishing of fj^ and -^rgf^ dX r = a and r = R. This enables (3) to be written more 
explicitly as 



G ■' (x, X , w) = 2^ ^^ — (4) 



hnnX In 



with A-i^n i-^) normalized according to the prescription 



d^M^Li^) ■ MVi*n' W = Sx.X'Sij'S„,.^m'Sn^n' (5) 



which in turn is seen to imply the conditions 



and correspondingly for g^ and g^_^. It is to be noted that (4) applies both to the interior and exterior domains 

provided only that the appropriate eigenfrequencies uj^^ and eigcnfunctions Aj^^ are applied in each case. 

Having determined the relevant propagators the problem of computing the Casimir energy Ec can now be addressed. 
Using the energy density given in (1), the propagators (4), and the normalization conditions (5), the appropriate 
derivatives can be taken together with the limits x' ^ x and t' -^ t. This readily yields the result 






InX 



where the summation is to be taken over both the interior and exterior eigenmodes of the system. This sum has been 
carefully evaluated in ref.9 using an exponential cutoff with results entirely consistent with previous calculations. 
Thus it only remains to be determined whether similar conclusions follow from the stress tensor method. 

To pursue this issue it is to be noted that the interpretation of r'^'(x) as the k component of the force per unit 
area normal to a surface ai comes from the k component of the conservation law 



d^Tf"" ^ 0. 



Thus 



— [ dx T*^" 



dcJi T 



kl 



where the integration on the rhs is to be taken over the surface E which bounds the volume V . Since the Ihs is the 
time rate of change of the k component of the momentum, it can then be claimed that T'^'(x) is the correct force per 
unit area so long as the calculation is performed in cartesian coordinates. However, the claim has been made [7,10-12] 
that it is also valid in spherical coordinates, i.e., that the force per unit area F/A on the spherical surface is simply 
given by 



F/A= {0\[Trr{r = a - e) - r„(r = a + e)]|0). 



(6) 



This, of course, is an assertion which can be examined by direct calculation, a task to which attention is now directed. 
At r = a the vacuum expectation value of Trr can be written as {0\^{B^ — _E^)|0}. Application of the boundary 
conditions at r = a together with the identity 



Ei^'™(^' 



21 + 1 

47r 



yields the vacuum expectation value of T„-{r ~ a — e) as the angular independent form 



(1). 



dr 
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5r„K;» 



(2) 



(7) 



where cul^ refers to the internal eigenfrequencies. The corresponding result for {0\Trr{r = a + e)\0) is simply obtained 
from (6) by using the corresponding external eigenmodes together with the replacements /^ -^ /^ and g< -^ g>. 
In order to relate this result to the mode summation method, however, an explicit construction of the relevant 
eigcnfunctions is required. For the interior region r < a the normalization of the functions fj^ and g^ leads to 



fiii^, 



(1),) = 



2uj 
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In 



d^^'^"' 



(1) 



a) 
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and 



(2) 



5,<RtV) 



Upon using these results it is found that 



/(/ + !) 

, ,(2) 2 



(2). 
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^Ml™ 4(4^)^^ + 1)4 



(A) 



(8) 



ln\ 



with the sum to be performed over the interior eigenfrequencies. If a corresponding resuh were to obtain for r,.^ over 
the exterior of the sphere, it would then follow that 



F/A\a-e+F/A\a+,^-{ 



)Ec 






(9) 



in agreement with the claim of [7]. However, for the exterior of the sphere there is, of course, a contribution only from 
r = a while the normalization depends upon both a and _R, a circumstance which in fact prevents one from obtaining 
a result analogous to (8). To see this it is convenient to go to manifestly normalized forms by the replacements 



!T.M^'.r 



2uj 
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In 






fi:.i^l'Jr) 



and 
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Upon using these expressions in the evaluation of {0\Trr{r = a + e)|0) there results 



a Aira'' 
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where 
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and 
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Since the coefficients a^^' vanish in the large R limit, it follows that this implies a suppression of the contribution 
of the external modes in the calculation of the Casimir energy,thereby contradicting Eq.(8). Thus the stress tensor 
approach does not constitute a valid approach to this problem for the case of spherical boundaries. Similar conclusions 
are readily obtained for Dirichlet and Neumann boundary conditions. 

Although the conclusions reached here have been based on a very specific calculation, they can be placed in a more 
general context by using techniques of covariant differentiation. To this end one takes the conservation law for T'^" 
in the form 



Ttr = 



d^T'"' + \ ^ |r"'^ + ( " |t^" = 0. (10) 

For spherical coordinates the nonvanishing independent ChristofFel symbols are 

1/r, 



r 9 j [r (j) 



u 



-r, i , , r" = —r^sin^t 



-sinOcosO, and < „ > = cot9. 



e 

This yields for the r component of (10) that 

r"^ smO r 

where use has been made of the tracelessness oiT^". Since (0|T''"|0} is time independent and (0|r'"''|0) — (0|T^'^|0) 
0, it follows that 

I drdndr{r^{Q\T"'\{))) ^ f r^drdn(0\T°"\0}, 

and consequently 

„2 j„/r,n^OO/ 



r^(0|r„(r)|0)|;j == / r'dr{0\T"^{r)\0) (11) 

for arbitrary ri and r2. This clearly shows that the correct inference to be drawn from (10) is not the force equation 
(6) but rather the fact that the difference between the quantities r^ {Q\Trr{r)\0) when evaluated on the two bounding 
spherical surfaces is proportional to the total energy in that region. 

It should also be noted that this result can also be inferred directly from the cartesian tensor result 

di{0\T^^\0) =0. 

which implies that 

Xkdi(0\T'^'\0) - dix,{0\Tki\0) - (0|T'='^|0) 
= 0. 

Upon using once again the tracelessness of the energy momentum tensor Eq.(ll) clearly follows. 

It has been noted here that many of the existing calculations of the Casimir effect for a sphere have not used correct 
boundary conditions on the underlying propagators. More serious is the fact that those which have used the stress 
tensor method have applied an approach which fails in its application to nonplanar boundaries. In concluding this 
work it needs to be pointed out that there exists yet a third failure in many of these calculations. To illustrate this 
point it is convenient to refer to Eq.(3.14) of ref.[7] which gives the Casimir energy as 

i v^/ , N f°° d(uja) 



2a ^' './__ 2tt 



(zjiY , (zjir , i^h['^y , i^h\'^)"' 



ire and the / 
bring this result to the usual integral over Bessel functions of imaginary argument it is necessary to perform a ninety 



where the ji terms are associated with the interior of the sphere and the h\ terms with the exterior. In order to 



degree rotation of the contour. However, for the case of the exterior mode part of Ec such a rotation fails since it 
does not take into account the existence of poles of the Hankel functions in the lower half plane. This is a noteworthy 
reminder of the remarks made earlier concerning the fact that the outgoing spherical wave condition for r > a cannot 
be derived for Casimir propagators. 

Having displayed some of the problems associated with the calculation of the Casimir energy for the case of spherical 
boundaries it is of interest to note that calculations which employ the stress tensor method claim to obtain the same 
result as that found by direct mode summation. This issue is dealt with in an appendix which shows that when 
appropriate attention is paid to the issue of contour rotation one does not in fact obtain the usual result. 

APPENDIX 

It has been stated in this work that the discontinuity in the stress tensor across a boundary cannot be used to 
calculate the Casimir pressure on that surface. Since, however, such calculations invariably claim to obtain the usual 
result, it is of considerable interest to display explicitly the flaw in such calculations. In order to avoid inessential 
complications one can choose to deal only with Trr{r = a + e) for the case of the TE modes with corresponding results 
for Trr{r = a ~ e) and for the TM modes then following immediately. This is, of course, equivalent to considering 
the scalar field case with Dirichlet boundary conditions. For this case one finds in analogy to [7] and [12] that the 
relevant Green's function is obtained by solving the equation 

(- V^ - u^) G(x, x'; uj) ^ (5(x - x') 
subject to Dirichlet boundary conditions at r = a, R. Upon writing 



G(x, x'; c.) - ^ G^r, r'; c.)rr(^, 0W"*(^', </-') 



1=0 

with y^™ denoting the usual spherical harmonics, it follows that 

where W[f^ , ff"] denotes the Wronskian and r> (r<) denotes the greater (lesser) of r and r' . In terms of the spherical 
Bessel functions ji and n; one readily finds that f^ {i-^r^ ) and ff" {ujr^ ) have the form 

fi^iujr) ^ ji{ujr)ni{uja) - ii{uja)ni{ujr) 

and 

ff'iujr) ^ ji{Lor)ni{ujR) - ji{ujR)ni{ujr). 

This allows the partial wave Green's function to be written in the form 

ni[Loa)]i{ULiR) - ni(u;R)ji[uja) 

Using the fact that the zeros of the Wronskian determine the eigenmodes of the system Gi{r, r';t — t') is found to be 
given by 

where the sum is taken over all eigenmodes of the system corresponding to eigenfrequencies LOin- This form of the 
propagator is recognizable as the Fourier transform of the scalar version of (3) with explicitly normalized eigcnfunc- 
tions. An immediate consequence is that the calculation of Casimir energy by direct mode summation can now be 
carried out using the approach of [9] to obtain the usual result. In addition the analysis following Eq.(9) can be 
repeated, again finding that there is no contribution to {Q\Trr{r — a + e)|0} in the large R limit. 

Although this would seem to establish the inapplicability of the stress tensor method, it is instructive to continue 
the analysis to determine exactly at what point the latter method fails. To this end one notes that the evaluation of 
(0|Trr(r = a + e)|0) requires that derivatives with respect to both r and r' be taken at r = r' = a and an integration be 
performed over w. Since the eigenmodes occur at real values of w, the path of integration must be carefully specified. 



As deduced from the boundary conditions (and also as stated clearly in [7]) the appropriate path is just above (below) 
the real axis for w > (w < 0). The usual expression for the Casimir force is then obtained if (a) the contour can 
be rotated ninety degrees counterclockwise and (b) if f{'{uja)) becomes proportional to the Hankel function h\ ' {uja) 

(2) 

(hi {uja)) in the upper (lower) half plane in the large R limit. Although the second condition (b) is satisfied, giving a 
R independent result which is formally identical to the usual expression for the Casimir energy, it is easy to see that 
the required contour rotation cannot be performed. To display this result one writes 

/OO - T 7 

which can be reduced to 

, , . ., > -« f°^ ■ u f'dii^R)^M^r) ~ ni(ujR)-fji(ujr) 

2a^ t=t'+e J_^ ji [LoRjni (wa) - ji [ujajni [ujR) 

While it is clear that the segment of the contour from lu ~ —oo to lu ^ can be rotated counterclockwise to the 
negative imaginary axis, no such rotation is possible for the segment lo = to lo ^ oo [16]. Conversely, in the case 
that the stress tensor is defined by the limit t = t' — e rather than t — t' + e, it becomes possible to rotate that part of 
the contour along the positive real axis to the positive imaginary axis. This choice, however, precludes the possibility 
of a legitimate rotation of the contour along the negative real axis. This completes the proof that the stress tensor 
does not in fact yield the usual expression for the Casimir stress on the sphere. The demonstration in the body of this 
paper was based on the method of expansion in eigenfunctions and in this appendix on the Wronskian formulation 
of the Green's function, the latter being the one more commonly used in published calculations of the stress tensor. 
The results of the two approaches are identical (note the remark [16]) and show that stress tensor calculations in 
curvilinear coordinates cannot be expected to yield the correct Casimir energy. 
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